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In our previous works the effective potential approach for composite operators has been generahzed 
to non-zero temperature in order to derive the analytical equation of state for pure 5'?7(3) Yang- 
fSJ , Mills fields without quark degrees of freedom. In the absence of external sources this is nothing but 

the vacuum energy density. The key element of this derivation is the introduction of a temperature 
■ dependence into the expression for the bag constant. The non-perturbative part of the analytical 

04 , equation of state does not depend on the coupling constant, but instead introduces a dependence on 

the mass gap. This is responsible for the large-scale dynamical structure of the QCD ground state. 
^ ■ The perturbative part of the analytical equation of state does depend on the QCD fine-structure 

coupling constant as. Here we develop the analytical formalism, incorporating the perturbative part 
in a self-consistent way. It makes it possible to calculate the PT contributions to the equation of 
state in terms of the convergent series in integer powers of a small as. We also explicitly derive and 
numerically calculate the first perturbative contribution of the Os-order to the non-perturbative 
part of the equation of state derived and calculated previously. The analytic equation of state 
f~ l ■ or, equivalently, the gluon pressure is exponentially suppressed at low temperatures, while at the 

temperature T — Tc = 266.5 MeV it has a maximum, if divided by It demonstrates a highly 

non-trivial dependence on the mass gap and the temperature near Tc and up to approximately 
(3 — 4)Tc. At very high temperatures its polynomial character is confirmed, containing the terms 
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proportional to T^ and T with a non-analytical dependence on the mass gap. 
PACS numbers: ll.lO.Wx, 12.38.Mh, 12.38.Lg, 12. 38. Aw 



INTRODUCTION 



D 
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] Up to now, lattice QCD remained the only practical method to investigate QCD at finite temperature and baryon 
density [H-Q- Recently it underwent a rapid progress f [cl-fToj and references therein). However, lattice QCD is 
primarily aimed at obtaining well-defined corresponding expressions in order to get realistic numbers for physical 
quantities. One may therefore get numbers and curves without understanding what is the physics behind them. Such 
an understanding can only come from the dynamical theory, which is continuous QCD. For example, any description 
of the quark-gluon plasma (QGP) 0, Q has to be formulated within the framework of a dynamical theory. The 
need for an analytical equation of state (EoS) remains, but of course it should be essentially non-perturbative (NP), 
] approaching the so-called Stefan-Boltzmann (SB) limit at very high temperatures Ih, J^]. Thus the approaches of 
analytic NP QCD and lattice QCD to finite-temperature QCD do not exclude each other; on the contrary, they should 
be complementary. This is especially true at low temperatures where the thermal QCD lattice calculations suffer from 
big uncertainties (see above-cited papers), while any analytic NP approach has to correctly reproduce thermal QCD 
lattice results at high temperatures. 

The formalism we use to generalize it to non-zero temperature is the effective potential approach for composite 
operators [l^. In the absence of external sources it is nothing but the vacuum energy density (VED). This approach 
is NP from the very beginning, since it deals with the expansion of the corresponding skeleton vacuum loop diagrams, 
and thus allows one to calculate the VED from first principles. The key element in this program was the extension 
of our paper [l^ to non-zero temperature [TEj . This makes it possible to introduce the temperature-dependent bag 
constant (pressure) as a function of the mass gap which is responsible for the large-scale structure of the QCD ground 
state. The confining dynamics in the gluon matter (CM) will therefore be nontrivially taken into account directly 
through the mass gap and via the temperature-dependent bag constant itself, but other NP effects will also be present. 
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Being NP, the effective approach for composite operators, nevertheless, makes it possible to incorporate the thermal 
perturbation theory (PT) expansion in a self-consistent way. 

The main purpose of this article is the derivation of the analytical EoS for the GM, i.e., a system consisting purely 
of SU{3) Yang-Mills (YM) fields without quark degrees of freedom. Its NP part which solely depends on the mass 
gap has been evaluated in [15]. Its PT part together with the mass gap depends on the QCD fine-structure constant 
as ■ Here we are going to formulate and develop the analytic formalism which makes it possible to determine the PT 
part of the YM EoS or, equivalently, the GM EoS in terms of the convergent series in integer powers of as. So this 
allows to calculate the PT contributions termwise in all orders of a small as ■ We also explicitly derive and numerically 
calculate the first PT contribution of the ag-order to the NP part of the GM EoS derived and calculated previously 
[l5| . The low- and high-temperature expansions for the GM EoS have analytically been evaluated. 

The present paper is organized as follows. In section II the general expressions for the gluon pressure as a function 
of temperature are present. In section III all results for the NP part of the gluon pressure are collected and briefly 
explained. Thus in the short sections II and III we describe the results obtained previously in [13, [iBl j and here we 
present them only for the reader's convenience in order to have a general picture in hand. In section IV the analytic 
thermal PT is developed, while in section V we prove that it makes it possible to find the PT part of the gluon 
pressure in terms of the convergent series in integer powers of a small as- The gluon pressure up to the a^-order 
term is derived in section VI. In section VII we present our numerical results, where our conclusions are also given. 
And finally, in appendices A and B the low- and high-temperature expansions for the gluon pressure are analytically 
evaluated, respectively. 

II. THE GLUON PRESSURE AT NON-ZERO TEMPERATURE 

In the imaginary-time formalism [l^ [l6l . [l7| all the four-dimensional integrals can be easily generalized to non-zero 
temperature T according to the prescription (let us remind that in the present investigation the signature is Euclidean 
from the very beginning) 



dgo 
(2^) 



T 



+ CXD 

E 



9 = q + 90 = q + =0J^ + i^i, = 2mTT. 



(2.1) 



In other words, each integral over qo of a loop momentum is to be replaced by the sum over the Matsubara frequencies 
labeled by n, which obviously assumes the replacement — >■ a;„ = 2m:T for bosons (gluons). 
Introducing the temperature dependence into the gluon pressure [l^, we obtain 



Pg{T) = Pnp{T) + Ppt{T) = Bym{T) + Pym{T) + Ppt{T), 



(2.2) 



where the corresponding terms in frequency-momentum space are: 



Bym{T) 



-10 
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dw T ^ 

n— — oo 



ln(l + 3a^^^(c.2,..2)) 
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Pym{T) = J / dw T V 

JO 



\n[l + \a^«^{u:^.:l) 



(2.4) 



Ppt{T) = duju^T y 



In 1 



4 -f 3Q;-'"-^(aj^,a;„)/ 4 



(2.5) 



In frequency-momentum space the intrinsically non-perturbative (INP) and the PT effective charges become 



a [q ) = a (q ,w„) = a (w ,w„) = 



.INP/2 , ,2 



INP/, ,2 , ,2^ 



a;2 -t- oj?. 



(2.6) 
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and 



a^^iq') = a^^ie,u.l) = a^^ico^u^l) = (2.7) 
respectively. It is also convenient to introduce the following notations: 



T-^=/3, u; = V<f, (2.8) 

where, evidently, in all the expressions is the square of the three-dimensional loop momentum, in complete agree- 
ment with the relations (2.1). 

In Eq. (2.6) is the Jaffe-Witten mass gap [l^, mentioned above, which is responsible for the large-scale structure 
of the QCD vacuum, and thus for its INP dynamics. Recently we have shown that confining effective charge (2.6), 
and hence its /3-function, is a result of the summation of the skeleton (i.e., NP) loop diagrams, contributing to the 
full gluon self-energy in the g2 ^ regime. This summation has been performed within the corresponding equations 
of motion (and references therein). In more detail (including the interpretation of Eq. (2.6) and the explanation 
of all notations) the derivation of the bag constant as a function of the mass gap and its generalization to non-zero 
temperature has been completed in fl3l and [isj . respectively. Let us only note that we omit the subscript "s" in 
a^^^iq^) which has been used in |15| . 

The PT effective charge a^^ {q^} (2.7) is the generalization to non-zero temperature of the renormalization group 
equation solution, the so-called sum of the main PT logarithms [Mil]- Here Af,^.^ = 0.09 GeV^ H is the asymptotic 
scale parameter for SU(3) YM fields, and b = (ll/47r) for these fields, while the strong fine-structure constant is 
as = as{mz) = 0.1184 [2J]. In Eq. (2.7) q"^ cannot go below Ay^, i.e., Af^j^^ < q^ < oo, which has already been 
symbolically shown in Eq. (2.5). In our previous works the expression (2.7) was denoted as a^^{q'^). However, 

here we prefer to denote it as in Eg. (2.7), leaving the notation a^^{q^) for the asymptotic freedom (AF) relation 
a^^{q'^) = l/b\ii{q'^ /Aym) itself [l9l - l22| . One can recover it from Eq. (2.7) in the q^ ^ oo limit. 

The NP pressure Pnp{T) = Bym{T) + Pym{T) and the PT pressure Ppt{T), and hence the gluon pressure Pg{T) 
(2.2), are normalized to zero when the interaction is formally switched off, i.e., letting = = 0. This means that 
the initial normalization condition of the free PT vacuum to zero [l3l - [l5j holds at non-zero temperature as well. 



III. Pnp{t) contribution 

One of the attractive features of the effective charge (2.6) is that it allows an exact summation over the Matsubara 
frequencies in the NP pressure Pnp( T) g iven by the sum of the integrals (2.3) and (2.4). Collecting all analytical 
results obtained in our previous work |15| , we can write 

6 ,2.,^. 16. 



and 



Pjvp(T) = Bym{T) + Pym(T) = — A^Pi(T) + -.T[P2{T) + P^{T) - Pi{T)], (3.1) 



^i(^) = / d--^^, (3.2) 



while 



•''-'off 

PsiT) = [ duj uj^ In (l - e-^'^' 



PilT) = / dw w2 In (1 - e~ 
Jo 



(3.3) 



where LOef / = 1 GeV and the mass gap A2 = 0.4564 GeV^ for SU{3) gauge theory are fixed [l5| . Then uj' and u 
are given by the relations 
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3A2 



VSA = 1.17 GeV, 



(3.4) 



and 



cj= ya;2 + ^A2 ^ ^oj^ + fnl^^, me// = = 0.585 GeV, (3.5) 

respectively. Pnp{T) is shown in Fig. 1. It is worth reminding that in the NP pressure (3.1) the bag pressure Bym{T) 
(2.3) is responsible for the formation of the massive gluonic excitations oj' (3.4), while the YM part Pym{T) (2.4) is 
responsible for the formation of the massive gluonic excitations (D (3.5). 
Concluding, let us note that the so-called gluon mean number [l2| 



Ng^Ng{p,c,) = -^^, (3 = T-\ (3.6) 

which appears in the integrals (3.2)-(3.3), describes the distribution and correlation of massless gluons in the GM. 
Replacing by u) and uj' we can consider the corresponding gluon mean numbers as describing the distribution and 
correlation of the corresponding massive gluonic excitations in the GM, see integrals PsiT) and P4{T) in Eq. (3.3). 
They are of NP dynamical origin, since their masses are due to the mass gap A^ . All different gluon mean numbers 
range continuously from zero to infinity. We have the two different massless excitations, propagating in accordance 
with the integral (3.2), conventionally denoted as loi, and the first of the integrals (3.3), conventionally denoted as 
UJ2- However, they are not free, since in the PT A^ limit they vanish (the composition [P2{T) + PsiT) — P4{T)] 
becomes zero in this case). So the NP pressure describes the four different gluonic excitations. The gluon mean 
numbers are closely related to the pressure; in particular, its exponential suppression in the T — ^ limit and the 
polynomial structure in the T — >■ oo limit is determined by the corresponding asymptotics of the gluon mean numbers. 
For the explicit evaluation of the low- and high-temperature expansions for the NP pressure (3.1) see appendices A 
and B, respectively. 



IV. ANALYTIC THERMAL PT 



Let us begin here reminding that we were able to perform the summation over the Matsubara frequencies analytically 
(i.e., exactly) for the NP part of the gluon pressure (2.2). To do the same for its PT part (2.5) is a formidable task. 
The only way to evaluate it is the numerical summation over n and the integration over uj, which is beyond our 
possibilities at present (if it is possible at all). Our primary goal in this article is threefold. Firstly, to develop 
the analytical formalism which makes it possible to calculate Ppt{T) (2.5) termwise in integer powers of a small as- 
Secondly, to calculate explicitly the PT contribution of the as-order. Thirdly, to derive the low- and high-temperature 
expansions for the gluon pressure. 

For the first goal, it is convenient to re-write the integral (2.5) as follows: 



Q /-OO +QO 

Ppt{T) = ^ / dujLo^T Y 



(4.1) 



where 



^ ' 4 + 3«^^^(w^c.2) iM{oj^,ujl)l + as\nzn' ^ ' 

with the help of the expressions (2.6) and (2.7). Here 



and (I>2 is given in Eq. (3.5). Let us also note that in these notations 
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1 + Q!c In 2„ 



(4.4) 



There is an interesting observation concerning the argument x{uj'^ ,ll>'^) of the logarithm ln[l + x(uj'^ ,lo'^)] in the 
integral (4.1). At its lower limit to — Aym and n — the argument (4.2) numerically becomes 



4A2,^^ + 3A2 



= 0.0185, 



(4.5) 



and the numerical values of A^, as and AyM given above have already been used. The argument of the logarithm is 
really small (it is an order of magnitude smaller than as itself), and it will become even smaller and smaller with uj 
and n going to infinity. This means that the logarithm ln[l + x{uj'^ ,0:'^)] in the integral (4.1) is legitimated to expand 
in powers of small x{ui'^,u!'^) at any n and in the whole range of the integration over uj, that's 00 > lu > Aym- Doing 
so, one obtains 1251 



\n[l + x{Lo\ial)] ^-^x"^{uj^^l), x{u^\^l) « 1. 



(4.6) 



Extracting the first term in the expansion (4.6), the integral (4.1) becomes 



Ppt{T) 



J A 



j.00 +00 

A V A r ^ — 



n— — 00 



o ^ ( 1 \ m 

r / 2 2\ / M y ^) mi 2 2\ 

:W„) - -a(z„)J - X (UJ ,UJJ 



m=2 



From now on it is instructive to separate the two terms in the integral (4.7) as follows: 



Ppt{T) - Ppt{A^;T) + P'pt{T), 



where 



Q -\-oo 



n— — 00 



A/(w2, w^) 1 + Q!s In z„ 



on account of the relations (4.2)-(4.4), and 



(4.7) 



(4.8) 



(4.9) 



8 f°° 



+00 



(-1)^ 



-x-^ico^^l) 



(4.10) 



A. PpT(A^;r) contribution 

The principal difference between these two terms (4.9) and (4.10) is that the former term vanishes in the formal 
PT A^ = hmit, while the latter one survives it. Let us consider the first term in more detail. The function 
(1 + as lnz„)~^ in the integral (4.9) can formally be replaced by the expansion in integer powers of as, namely 



(1 + as In z„ 



fe=0 



{-l)''a':\n' z„, 



(4.11) 
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which converges to the corresponding function in the whole range h-yu < uj < oo and at any n for ja^ Inznj < 1. The 
integral (4.9) can equivalently be re-written as follows: 



which makes it possible to present it as a sum in integer powers of a^, namely 



(4.12) 



PpT(A2;T) = ^a^-F,(A2;T), 

fc=i 



where 



(4.13) 



+ CXD 



For example, the first term Pi(A^; T) explicitly looks like 



1 



Af(w^tj2 



-(-1)^-1 ln'=-iz„ 



(4.14) 



Q r°° 1 



n— — oo 



2 tj2) 



(4.15) 



where M(a;^,a;^) = + oj^, and itself is given in the relation (3.5). In this integral the summation over the 
Matsubara frequencies can be performed analytically (i.e., exactly) with the help of the expression explicitly given in 
fisj l . Omitting all the derivation and dropping the /3-independent terms [l3| , one obtains 



27r2 Jav,. we/^^-l 



(4.16) 



B. P'pt{T) contribution 

On account of the relations (4.2)-(4.4), the integral (4.10) becomes 



P'pAT) 



./Ay 



+ 00 

n— — oo 



.■m=2 



(1 + as In ZnY' 



where 



(4.17) 



(-1)™3™ {uj^+ujI)"' 



m4™ M"^{uj'^,ujI)' 
In complete analogy with the expansion (4.11) one gets [1^ 



(4.18) 



(1 + as InZnY 



Ck{m)a^s 111'' 



where 



(4.19) 



co(m) = 1, Cp{m) 



■^(fcm-p + fc)(-l)'=Cp_fc, p>l. 

k=l 



(4.20) 
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What has been said in connection with the expansion (4.11) is valid for the expansion (4.19) as well. So on its account, 
the integral (4.17) can equivalently be re- written as follows: 



6„(w2^u;^)a™^cfc(m)a^ln'=2„ 



.m=2 



fe=0 



(4.21) 



Let us consider the coefficients 6to(w^,w^) (4.18) in more detail. Noting that 



(u;2+a;^) = M(u;^c^^)- 



these coefficients can be present as follows:. 



&m(w ,w„) = - - 



3\"' 1 {io^+ulr 



a) mM"(w2,w2) 



6\ 1 X ^ / m 



4 / m 
4/ m 



m — 1 



and 



m\ m(TO — l)...(m — p + 1) /m 





m 



= 1 



are the binomials coefficients. Substituting the expression (4.23) into Eq. (4.21), it becomes 



(4.22) 



(4.23) 



(4.24) 



P'p^iT) = Ppt{T) + P'priA-'-^T), 



(4.25) 



where 



oo 

E 

.m=2 



4/ m 



^^Cfe(m)ajln'= 2;„ 



(4.26) 



and 



o /■oo 



+ 00 



.m=2 



fe=0 



(4.27) 



with 



It is convenient to present the integral (4.26) in the following way 



(4.28) 



Q i^? 



^ V 4/ m + 2 

.m=0 ^ ^ fe=0 



^Cfe(m + 2)a^ln'=^„ 



(4.29) 
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which shows exphcitly that it is an a^-order term. 

It is also convenient to present the integral (4.27) in the same way, namely 



duj to T 



+00 

E 

n— — oc 



E 

_m=0 



a" 



k=0 



00 

, 2\ \ " 



rn 



2)a^ In'^ 



(4.30) 



with 



m+l / I r)\ / \ 

p=o \ " / V 

r, \ 2 m+l / I 9\ / ^ \ ™~P 



p=0 

2 



3 /v^ , 

' '"m+2V 

(4.31) 

Then the previous integral (4.30) becomes 



P^^(A2;T) = 



9a^A^ 
V^47r 



+ 00 



00 
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m— 



(-3/4) 



-p: 



00 



m m 

-'(")(A2)> cfc(™ + 2)ann«z„ 



fc=0 



(4.32) 



and in the formal PT A^ = limit it is zero as the whole expansion (4.13). The integral (4.32) is really of the order 
of a^. However, this order numerically is much smaller than the corresponding order term in the expansion (4.13). 
This will be true for any corresponding orders in the expansions (4.13) and (4.32) because of the initial condition 



x(cj , ( 
other. 



^ 1 in Eq. (4.10). Obviously, the structures of all expansions (4.13), (4.29) and (4.32) differ from each 



V. CONVERGENCE OF THE PERTURBATIVE SERIES 



The convergence of the power series (4.10) over m is obvious, since it comes from the expansion of the logarithm 
(4.6). In turn, this means that the series in integer powers of as (4.29) and (4.32), due to the relation (4.25), 
are also convergent. At the same time, the convergence of the series (4.13) over k is not so obvious and requires 
some additional investigation. The series (4.13) can be formally considered as a power series over small as with the 
coefficients Pfc(A^;T). The QCD fine-structure coupling constant correctly calculated at any scale is indeed always 
small. Let us remind that in this paper we use its value calculated at the Z-boson mass, namely as — as{mz) — 0.1184 
[23 ]. For the convergence of the power series (4.13) it is necessary to show that its radius of convergence r is bigger 
than any possible value of as ■ The radius of convergence of the power series can be calculated in accordance with the 
Cauchy-Hadamard (CH) theorem as follows: 



lim 



Pfe(A2;T) 

if this limit exists. Substituting the corresponding expressions from Eq. (4.14), one obtains 



(5.1) 



lim 



r 

JA 



duj 



+00 



n— 



1 In*:-! 7 



1. 



(5.2) 



Thus the series (4.13) converges absolutely for as < r — 1 and converges uniformly on every compact subset of 
{tts .as < r}. Roughly speaking this means that any further calculated term in integer powers of as will be smaller 
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than the previous one (at least by one order of magnitude). Let us also note that now there is no need in the restriction 
\as lnz„| < 1, only < 1 is required, which always holds. 

Summing up all the PT contributions, the PT pressure (2.5) finally becomes 



PpTiT) = PpTiA^;T) + Pi^TiT) = [PpTiA^:T) + P'pTiA^;T)] + Ppt{T). (5.3) 

Due to this decomposition a few things should be made perfectly clear. All three integrals (2.3), (2.4) and (2.5), 
contributing to the gluon pressure (2.2) are the corresponding parts of the initial skeleton loop integral derived in 
[Til l. The first two integrals depend on the mass gap only, so they are truly NP, indeed. We call the integral 
(2.5) or, equivalently, (5.3) the PT one because it explicitly depends on as, but this is only a convention. Its two 
terms Ppt{A'^;T) (4.13) and P'prp{A^;T) (4.32) depend on the mass gap as well (they vanish in the PT A^ = 
limit). Its third term Ppt{T) (4.29) begins with the a^-order correction, and does not depend on the the mass gap 
at all. However, it is an infinite sum of pure PT contributions, i.e., in fact this is the corresponding part of the 
skeleton loop term (5.3). The terms Ppt{A^]T) and Ppj,{A^]T), being by themselves the two other corresponding 
parts of the skeleton loop term (5.3), can be considered as the ag- and a^-order corrections to the skeleton loop term 
Pym{T) (2.4), since they depend on the massive gluonic excitation w. Let us emphasize that there cannot be any 
PT corrections to the bag constant Bym{T), since its skeleton part (2.3) has been defined from the very beginning 
by the subtraction of all the types of the PT contributions [3, [iBl- AH three expansions (4.13), (4.29) and (4.32) 
analytically depend on as, and they can be calculated termwise in integer powers of as- The convergence of the series 
in integer powers of a small as derived here has been confirmed. It guarantees that any further PT contribution in 
powers of as will be numerically smaller than the previous one. It is worth emphasizing that in any case none of 
the PT contributions, and hence none of their sum, can be numerically bigger than the SB term, which describes the 
thermodynamic structure of the GM at high temperatures, see Figs. 1 and 7 below. 

Here a few general remarks are in order. In the initial thermal PT QCD the dependence on as is non-analytical. 



i.e., the expansion contains its fractional powers, a^lnas, etc. (see, for example [5, 12?, |17|, |26|, |27[ and references 
therein). This leads to the di verg ent series in the thermal PT QCD, though each term up to the Inas-order has 
been calculated correctly. In [l^ it has briefly been explained why this effect occurs there. This problem does not 
occur in our formalism, which is the NP from the very beginning as underlined above. Apparently, the reason is that 
it allows to calculate the Ofg-dependent corrections, which by themselves are the corresponding skeleton parts of the 
skeleton loop integral (2.5), to the NP term (2.4) which is the skeleton part of the initial skeleton loop integral (2.2). 
In other words, the re-summation of an infinite number of the corresponding contributions has been already done for 
all of them. So our approach allows to develop the PT series in integer powers of as for the NP quantities, i.e., the 
coefficients of these infinite series are the NP quantities. Let us emphasize that these conventionally called PT series 
are, in fact, the so-called cluster expansions [18|, which are convergent series and depend on the coupling constant 
not in elementary way, like all our expansions discussed here. This is completely different from the thermal PT QCD, 
which in the best case is of the asymptotic-type expansion. It is dealing with the pure PT contributions from the very 
beginning, and the necessary re-summation of the hard thermal loops (HTL) leads finally to the non-analytical 
dependence on as- At the same time, the formalism developed here makes it possible to calculate the as-dependent 
contributions to the gluon pressure in terms of the convergent series in integer powers of a small as. This seems to 
resolve the above-mentioned long-term problem. 

However, instead it may produce another problem, namely the question of double-counting in as- In the present 
investigation there is no double-counting in as, since all the convergent PT series (4.13), (4.29) and (4.32) are of 
different structure, as emphasized above. The series (4.29) and (4.32) begin with the a^-orders, and therefore their 
numerical contributions are very small, since the series are convergent. The convergent series (4.13) begins with the 
as-order, which is the only one to be numerically calculated below. The double-counting problem will indeed arise 
when we will include the SB term into the full EoS, since it should be approached in the AF way, which involves all 
powers of as . So the question how far do the pure PT corrections overlap with the NP ones cannot be answered here. 
But the problem is fixed, and we will address it and clarify the situation in the forthcoming paper. 



VI. THE GLUON PRESSURE Pg{T) 
Summing up all contributions, the gluon pressure (2.2) thus finally becomes 

Pg{T) = PNpiT) + Ppt{T) = Pnp{T) + [Ppt(A2; T) + P'p^{A^;T)] + Ppt(T). (6.1) 

In general, both expansions Ppt{A^;T) (4.13) and Pprp{A^;T) (4.32) are to be considered as producing the corre- 
sponding PT corrections to the leading NP part Pnp{T) (3.1) of the GM EoS (6.1). At the same time, the pure 
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FIG. 1: The gluon pressure (6.3), the NP pressure (3.1) and the as-dependent PT pressure (6.4), all scaled (i.e., divided) 

by r''/3, arc shown as functions of T/Tc- Effectively, all curves have maxima at Tc = 266.5 MeV (vertical solid line). The 
horizontal dashed line is the SB constant 3P<;s(T)/T*' = (24/45)7r^ 5.26. 



PT term Ppt{T) (4.29) is to be considered as producing the PT corrections to the leading PT contribution which is 
nothing but the above-mentioned SB term. However, due to the normalization condition of the free PT vacuum to 
zero (i.e., normalization of Pg{T) to zero when the interaction is switched off), it is not explicitly present in Eq. (6.1) 
(see discussion below) . So none of pure PT corrections has to be calculated unless the leading SB term is 
restored to Eq. (6.1) in a self-consistent way. In this connection, let us note that the SB term or, equivalently, 
the pressure of a gas of massless free gluons (ideal gas) can be considered as the = 1-ordcr pure PT contribution to 
the full pressure. That is why we start the numerical calculation of the PT contributions to the gluon pressure (6.1) 
from its first non-trivial order, namely the ag-order in the expansion (4.13). As emphasized above, it is the a^-order 
contribution to Pym{T) in the NP term (3.1), which is already present in the GM EoS (6.1). Then it looks like 

Pg{T) = Pr,p{T) + P^^(T) + 0(a2), (6.2) 

where PpriT) = a,,Pi(A2;r) and for Pi(A'^;T) see Eq. (4.16). Omitting the terms of the 0(a2)-order, for con- 
venience, it is instructive to explicitly gather all our results from the relations (3.1)-(3.3) and (4.15) for the gluon 
pressure (6.2) once more as follows: 

fi 1 6 

PgiT) = — A2Pi(T) + —T[P2{T) + PsiT) - P^T)] + P^^^T), (6.3) 
where the integral Ppj,[T) is 

PMT) = a.x-A^y^^^d..^-^, (6.4) 

while all other integrals Pn{T), n= 1, 2, 3, 4 are given in Eqs. (3.2) and (3.3). This form is convenient for the numerical 
calculations. Let us note that when the interaction is formally switched off, i.e., letting = = 0, the composition 
[P2{T) + PsiT) - P4{T)] becomes identical zero (see Eqs. (3.3)) and thus Pg{T) itself. The gluon pressure (6.3) and 
its first PT contribution of the a^-order Ppj.(T) (6.4) are also shown in Fig. 1. However, it is worth emphasizing once 
more that, in fact, the term P|,y(T) is NP, depending on the mass gap A^, which is only suppressed by the a^-order. 

VII. RESULTS, DISCUSSION AND CONCLUSIONS 
A. Results 

From our numerical results it follows that all three pressures Pg{T), P]\!p{T) and Pprp{T) as functions of T/T^ 
effectively have maxima at the same "characteristic" temperature T = = 266.5 MeV, when they are scaled (i.e., 
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FIG. 2: The gluon pressure (6.3), the NP pressure (3.1) and the as-dependent PT pressure (6.4), all properly scaled in GeV 
units, are shown as functions of T/Tc- 
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FIG. 3: The gluon pressure (6.3), the NP pressure (3.1) and the as-dependent PT pressure (6.4), all properly scaled in GeV^ 
units, are shown as functions of T/Tc. 
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FIG. 4: The gluon pressure (6.3), the NP pressure (3.1) and the as-dependent PT pressure (6.4), all properly scaled in GeV^ 
units, are shown as functions of T/Tc. 
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FIG. 5: The composition (B25) scaled by is shown as a function of T /Tc- It approaches zero from below. This means that 
it does not contain the constant SB term, i.e., the SB terms have already been exactly canceled far below STc. 
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FIG. 6: The high temperature asymptotics of the gluon pressure (6.3), the NP pressure (3.1) and the Q^-dependent PT pressure 
(6.4) in GeV^ units are shown as functions of T/Tc. At T ^ 23Tc the NP pressure Pnp{T) goes below the PT pressure Ppt{T). 



divided) by T*/3, see Fig. 1. From this figure it clearly follovirs that the gluon pressure (6.3) will never reach the SB 
constant (24/45)7r^ — 3Psb{T)/T'^ limit at high temperatures. That is not a surprise, since the SB term has been 
canceled in the gluon pressure from the very beginning due to the normalization condition of the free PT vacuum to 
zero (13I4T5I (see discussion below). On the other hand, from Fig. 1 it follows that all excitations of the NP dynamical 
origin, described by Pmp{T) and Ppj^iT), and hence by Pg(T) itself, do not survive in the very high temperature 
T — 7> 00 limit. From Fig. 1 it is also clear that the first PT correction (6.4) is indeed smaller than the NP term 
(3.1) in the moderately high temperature range up to approximately (3 — 4)Tc 15]. At the maximum the PT term 
numerically is 0.236, while the NP term is 2.63, i.e., the former is by one order of magnitude smaller than the latter 
one. Due to the chosen scaling both pressures rather rapidly approach zero in the limit of very high temperatures. 

It is worth noting here that the PT term (6.4) of the gluon pressure (6.3) describes the same massive gluonic 
excitations ui (3.5) which have already been created by the NP dynamics in the GM. However, they are suppressed 
by the a^-order, and therefore can be considered as new massive gluonic excitations in the GM. We will denote it as 
as • ui. So in total the gluon pressure (6.3) describes the five different massless and massive gluonic excitations within 
our approach. In any case, the different types of the massive and massless gluonic excitations of the dynamical origin 
will necessary appear at non-zero temperature [1, [HI, .28 30] (and references therein). 

Let us briefly discuss the asymptotic properties of all three pressures in more detail. Below Tc all pressures are 
exponentially suppressed in the T — >■ limit, see Fig. 1. This is explicitly shown analytically in appendix A by 
considering the integrals (3.2)-(3.3) and (6.4) in the above-mentioned limit. Of course, this suppression is related 
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FIG. 7: The high temperature asymptotics of the gluon pressure (6.3), the NP pressure (3.1) and the -dependent PT 
pressure (6.4) in GeV' are shown as functions of T/Tc. At T = 23Tc the NP pressure Pnp goes below Ppt{T). The SB 
pressure Psb{T) = (8/45)7r^r'', formally extended up to zero temperature, is also shown. 



to the low-temperature asymptotics of the gluon mean number (3.6), as mentioned above. The high-temperature 
(T — >■ oo) expansion is explicitly evaluated in appendix B. At moderately high temperatures up to approximately 
(3 — 4)Tc the exact functional dependence on T remains rather complicated. It cannot be determined by the analytical 
evaluation of the above-mentioned integrals - only numerically as shown in Fig. 1. This non-trivial T-dependence 
can also be seen in Figs. 2, 3 and 4. In each of these figures all three pressures are scaled in the same way. Fig. 5 can 
be interpreted as clear diagrammatic evidence of the exact cancelation of the SB terms close to T^. It is analytically 
shown in Eq. (B25), i.e., it occurs within the composition IQTMiiT) /ir'^T'^ . In the NP pressure (3.1) the exact 
cancelation of the mass gap terms I^^T^ also occurs somewhere rather close to Tc, see Fig. 3. For the analytical 
evaluation of this phenomenon see appendix B, in general, and the high-temperature expansion (B24), in particular. 
As a result, the NP pressure (3.1) will scale as T, while the PT pressure (6.4), and hence the gluon pressure (6.3), will 
continue to scale up to the leading order as ^ a^A^T^ in Eqs. (B33) and (B34), respectively. Thus both pressures 
will approach the same constant in the limit of high temperature in Fig. 6, but very slowly since the term ^ is 
suppressed by the as order. At T = 23rc the NP pressure goes below the PT one, see Figs. 6 and 7. In general, 
all pressures are polynomials in integer powers of T up to at very high temperatures. The term ~ has been 
first introduced in the phenomenological EoS [3l| (see also (32l - [36j and references therein). On the contrary, in our 
approach both terms ~ and ^ T have not been introduced by hand. They naturally appear as a result of the 
explicit presence of the mass gap from the very beginning in the NP analytical EoS [l^ . 

B. Conclusions 

Our final conclusions are as follows: 

(i) . The effective potential for composite operators [13] may provide a new general analytic approach to QCD at 
non-zero temperature and density [l3 . [l5| . 

(ii) . It is essentially NP by origin, but may incorporate the thermal PT expansion. 

(iii) . We have developed the analytic thermal PT in the form of the convergent series, which made it possible to 
calculate the PT part of the gluon pressure termwise in integer powers of a small as ■ 

(iv) . We have shown that the PT contribution of the a^-order is numerically much smaller than the NP term in 
the range up to 23Tc, see Figs. 1-4 and 6-7. 

(v) . In the gluon pressure (6.1) the higher order terms in integer powers of a small as, which are determined by 
the convergent expansions (4.29) and (4.32), can be neglected. 

(vi) . All three pressures scaled (i.e., divided) by T^/Z have maxima at the temperature T = Tc = 266.5 MeV. Their 
low- (below Tc) and high-temperature (above Tc) expansions have been evaluated in appendices A and B, respectively. 

(vii) . In the low-temperature T — >■ limit all three pressures are exponentially suppressed (Fig. 1) due to the 
corresponding asymptotic of the gluon mean number. In the low-temperature expansion of the NP part (see Eq. A28)) 
a non-analytical dependence on the mass gap appears in terms ~ (A^)^/^r^ ~ AT^. 
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(viii) . The complicated mass gap- and T-dependence of all three pressures near and up to approximately 
(3 — 4)Tc is seen in Figs. 1, 2, 3, 4. In particular, the behavior of each pressure is not exactly power-type. This 
especially clearly follows from Figs. 2, 3 and 4. 

(ix) . The polynomial character of the high-temperature expansions for all three pressures is confirmed due to the 
corresponding asymptotic of the gluon mean number. 

(a) . For the NP pressure (3.1) it contains only terms ^ T, and some of them 
depend non-analytically on the mass gap, namely ~ (A^)^/^T ~ A^T. 

(b) . For the a^-dependent FT contribution (6.4) it contains the terms ~ 
and ~ T with a non- analytical dependence on the mass gap as above. 

(c) . For the gluon pressure (6.3) it contains both types of terms. 

(x) . In low- and high-temperature expansions, a non-analytical dependence on the mass gap occurs, as described 
above, but it is not an expansion parameter like is a^. 

(xi) . In the gluon pressure (B34) the mass gap term ^ A^T^ is explicitly present in the whole temperature range, 
though it is suppressed by the as-order. 

(xii) . The FT part (B26) dominates over its NP counterpart (Bl) in the limit of very high temperature starting 
from T = 23Tc, see Figs. 6 and 7. It is expected from the general point of view. This underlines once more the 
importance of the calculated here the as-dependent PT pressure. 

(xiii) . Our analytical derivations in appendices A and B are in complete agreement with our numerical results 
shown in Figs. 1-7 and vice-versa. 

(xiv) . In addition to the two massless wi, uj2 and the two massive w', uj gluonic excitations described by Pnp{T) 
(3.1), we have also one new massive excitation as • u) described by Ppj.(T) (6.4). All these excitations are of the NP, 
dynamical origin, since in the PT A^ — limit they vanish from the spectrum. 

The NP part (3.1) determines the thermodynamic structure of the GM at low temperatures within our approach. 
It is uniquely fixed and numerically it is half of the SB value at (Fig. 1 and [IBl)- The main problem which remains 
to solve is how to include the massless free gluons contribution to the gluon pressure (6.3), in order to reach the SB 
limit at high temperatures in the GM. Above Tc it will be always much bigger than any other correctly calculated 
PT contributions to the gluon pressure, see Figs. 1 and 7. The above-mentioned problem is not an easy task due 
to the normalization condition of the free PT vacuum to zero. The SB term cannot simply be added to the gluon 
pressure even multiplied by the corresponding Q{{T/Tc) — l)-function. This will lead to the finite jump at T = Tc in 
the full pressure which is not acceptable. The full pressure is always a continuous growing function of temperature 
at any point of its domain. So its inclusion into the YM EoS (which describes the thermodynamic structure of the 
GM in the whole temperature range) should be done in a more sophisticated way (work is in progress). Only after 
completing this program we will be able to analyze and compare our results for SU{3) GM EoS with the lattice 
results in ,3Z - ,34 . ,37 39] and analytical ones in (_5,, 27, 29- 3i 35., 3M (^'^^ references therein). The comparison and 
agreement with the thermal lattice calculations is especially important, since it will make it possible to understand 
what is the physics behind them. Our ultimate goal is the description of the dynamical structure of the QGP within 
NP analytical approach to QCD at non-zero temperature and density developed and advocated here and in (l^. [l5|. 

However, let us underline in advance that the gluon pressure Pg{T) analytically investigated and numerically 
calculated in the present investigation will be the most important part of the full pressure. Just it will be mainly 
responsible for its NP physics and will provide the agreement with lattice data in the whole temperature range and 
especially near Tc- The above-mentioned addition of the SB term will only ensure the correct high temperature limit 
of the full pressure. And finally, we do not expect that the inclusion of quark degrees of freedom will somehow affect 
the convergence of the corresponding PT series. But this has to be investigated and checked after SU (3) YM EoS is 
fixed, i.e., when the above-mentioned program is completed. 
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Appendix A: Low-temperature expansion 

Let us be gin with noting in advance that ah exactly calculated integrals, discussed in appendices A and B, can 
be found in [25l lioj . This is also true for the asymptotics in the low- and high-temperature limits of those integrals 
which cannot be analytically derived. In order to evaluate a low-temperature expansion for the gluon pressure 

Pg{T)^PNp{T) + P^T{T), (Al) 
it is convenient to present the NP pressure as in Eq. (3.1), namely 

Pnp{T) = -^A2pi(T) + ^TM{T), (A2) 
where the integral Pi{T) (3.2), on account of Eq. (3.6), is 

Pi(r) = / dujuj N,{p,Lo) = / d..^^, (A3) 

and 

MiT)^P2{T) + P3{T)-PiiT), (A4) 

with 



P. 



21 



(T) ^ / duj uj^\n{l-e-'^'^) 
P3{T) = J dw a;2 In (l - e-^"^" 

poo 

Pi{T) ^ / dtJ a;2 In (1 - e"^'^) . (A5) 



In all the above-displayed integrals the variable y — e^^^ is always small, and hence y^^ — e^^ is always big, in 
the low-temperature limit T — > (/3 = T^^ oo). This is true for the exponents e~^'^ and e~^" as well. Then the 
gluon mean number Ng{f3,uj) in the integral (A3) can be approximated as Ng{f3,uj) ~ exp(— /3aj) up to the leading 
order. So Pi{T) becomes 

/•CXD roo 

Pi(T)==/ duj uj Ng{(3,uj) dujuie-^'^. (A6) 

The almost trivial integration yields 

Pi(r) ^ (r2 + ^^^^r)e-^, r^o. (A7) 

The integral P2{T) can be considered in the same way. Up to the leading order it becomes 

.2 „-8ui 



duj uj^ e-P'^, P^oo, (A8) 



pOO pC 

P2(r)=/ doj uj^\n{l--e-^") ^ ~ 
and integrating it, one obtains 

P2{T)^^[2T^ + 2Lu,ffT^+LulffT)e-^, T ^ 0. (A9) 
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The integral PsiT) up to the leading order looks like 



PziT) = j dio io^\n(l-e~^^j ^ - J 



iff 



oo, 



(AlO) 



and replacing the variable co by the variable w' in accordance with the relation (3.4), this integral becomes 



P3{T) 



— ?► oo, 



where 



e// - V y^eff 



{ujlff+a?), a = \/3A. 



(All) 



(A12) 



Unfortunately, this asymptotical expression (All) cannot be directly evaluated, like it has been done for the integrals 
(A6) and (A8). However, noting that the variable x = /uj''^ < 1, we can formally expand 



k=2 



(A13) 



Then from the integral (All) one obtains 



duj' uj''^ e-^^' + - A2 / dJ e-l"^' + Pf') (T), /? ■ 



(A14) 



where 



Pt\T) 



k=2 



1/2 
k 



(A15) 



Let us consider the last integral (A15) in more detail. Since the series over k are convergent in the interval of 
integration and the functions depending on k are integrable in this interval, these series may be integrated termwise 
[25(] . that is. 



Pt\T) 



E 

k=2 



1/2 
k 



iff e 
diu' 



V3A (w') 



(A16) 



Integrating it, one obtains 



k=2 



'^'iff 
\/3A 



(A17) 



and 



N^\t,uj') denotes the result of the integration over lo' in Eq. (A16) on the interval [-\/3A, ^ f], while the 

Ak 



function Nf' (T, uj') itself is 



m-l( l\va+2-2k 



(-/3) 



2/0-3 



(2fc - 3)(2A: - A). ..{2k - 2 - m) 
The series for the exponential integral function Ei(— /3cj') is [1^ 



{2k - 3) 



-Ei(-/3w')- 



(A18) 
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/ " (1 — ^V 



■ oo. 



(A19) 



If one chooses n = 2fc — 3 in the previous equation (i.e., neglecting the relative error of the approximation (A19)), it 
is easy to show that both terms in Eq. (A18) for N^''\t,ijj') cancel each other termwise for any k>2, and thus 

(A20) 



ivi'=^(T,a;')=0, fc = 2,3,4,... 



or, equivalently, 



Pf)(r) = 0, fc = 2,3,4,... 



(A21) 



-'3 I — n — ... . y^^-^j 

Going back to Eq. (A14) and easily integrating the first two terms, and taking into account the previous result, one 
comes to the following expansion 



-A^Te — + 
2 2 

The integral Pi{T) up to the leading order looks like 



P4(r) 



du} up" In (l — e' 



(A22) 



(A23) 



JO JO 

and replacing the variable lj by the variable ui in accordance with the relation (4.5), this integral becomes 

/.oo 

P4(T) ~ - / dw ^^(w^ _ (a/2)2) e-f^'^, /3 ^ oo. (A24) 

J{a/2) 

Noting again that the variable z = a^/Aco'^ < 1, we can formally expand 



v/(a)2 - (a/2)2) = u> 
Then from the integral (A24) one obtains 



l-(aV8^^) + f:('f)(-z)* 



(A25) 



F4(T) 



dQ oj^ e-^^ + / duj e-^^ + Pi''(r), /? ^ oo. 



(A26) 



/\/3A/2 8 Jv3A/2 

where due to the same formalism which has been used previously in order to get the result (A21), one can conclude 
that P^\t) = as well. Easily integrating the first two terms, one comes to the following expansion 



Pi{T) {2T^ + VSAT^ + -A^T)e-'^ + ^A^Te"^, T ^ 0. 

4 8 



Substituting all these expansions into the Eq. (A2), one obtains 



(A27) 



Pmp{T) ^ ^A2(r2+^,,,r)e-^-i|r[2r3 + 2..e//T^+c.,2,,T]e-^ 



+ 



16. 

24 9 ,9 



✓3A 



{2T^ + 2u'^ffT^ + uj'^ffT)e-^ - {2T^ + 2^/1 AT^ + 3A^T)e-^ 



eff V3& 

e T — e T 



16^ 



2T^ + V3AT^ + -A^T 



V3i\ 

e"— , T-)- 0. 



(A28) 
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Evidently, this is nothing but a low-temperature expansion for the NP pressure Pnp{T). Let us note that it contains 
a non-analytical dependence on the mass gap squared in terms ^ (A^)^/^T^ ~ AT^, but the mass gap is not an 
expansion parameter like a^. 

It is instructive to re-write this expansion as follows: 



Pnp{T) ^ FNPiA,T) + bi 



16 



"iff "cff _V3& %/3A 



PsBiT) 



16 



(A29) 



since for YM fields (see, for example [Tl|, [ij] and references therein) 



Psb{T) - i-n^T\ (A30) 
45 

Fnp{A,T) denotes the sum of all the terms which directly depend on the mass gap, so that Fnp{T,A = 0) = 
and bi = (ISO/tt*). So the propagation of massless gluons below Tc can be described by the SB-type term. However, 
it is exponentially suppressed in the T — > limit, as it should be. At T Tc its contribution can be numerically 
comparable with other contributions in Eq. (A28). That's no surprise that the massless gluons may be present in the 
GM at any temperature. Moreover, let us note in advance that the propagation of the massless free gluons below Tc 
(if any) cannot be described by the SB term itself. It should also be exponentially suppressed in the same way as it 
is shown in Eq. (A29). But it has to survive in the PT = limit, while the contribution (A29) vanishes in this 
limit (Wg^j = We// m this case, see Eq. (A12)). The SB term can describe the propagation of the massless free gluons 
only in the high temperatures limit above Tc (see appendix B below). 

Let us now consider Eq. (6.4), which in the T — > (/3 = T^^ — > oo) limit up to the leading order becomes 

PpAT) - I^A^ r du i^3-L^ ^ I^A^ r dw w^ le-'^-, /? ^ oo, (A31) 
and ui is given by the relation (3.5). Replacing the variable oj by the variable w, as in Eq. (A23), one obtains 



PpAT) ^ -iA' / dw ^{co^ ~ {am e-^^ /3 ^ oo, (A32) 

where We// — \/ ^ym ('^Z^)^! ^nd for a see Eq. (A12). Noting that the variable z = a^/Auj'^ < 1 in this case, we 
can use the expansion (A25) in order to obtain 



27r^ 



da; 



-A^ 



do; 



■ff 



OJ 



P^'\T) 



/3 oo. 



(A33) 



Due to the same formalism which has been used previously in order to get the result (A21), one can conclude that 
p{k) j-y-j = as well. Easily integrating the first two terms, one comes to the following expansion 



PpAT) - 



(r2 + rwe//)e-^ + ^A2Ei(-^^) 



(A34) 



Here and below the corresponding exponential integral functions are defined by Eq. (A19). Summing up the expansions 
(A28) and (A34), one obtains a low-temperature expansion for the gluon pressure (Al) as follows: 
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16 
+ -T 

27r2 



(2T3 



':ffT)e- 



(2T3 + 2V3Ar2 + 3A2T)e-^ 



e — e 



2T^ + \/3AT2 + -A^T 



(r2 + T^e//)e-^ + ^A2Ei(-!i:^) 



8 

T-i^ 0. 



e 



(A35) 



Let us note that the low- temperature expansion (A35) depends mainly on the effective massive "excitation" uj'^^^ = 

.^Ju^^J^rcfl, and does not depend on the effective massive "excitation" uj^ff = \J^tff + at all, while the 

dependence on the effective massive "excitation" Cjeff = \/h.YM + (''■*/2)^ is suppressed by a^. 

It is instriictive to use in the exponents of the previous expansion the following obvious relations: We// = 
i^iTc, ijj'^jf = v^Tc, a/SA = V2,Tc, 1^4 = (1/2)^3, Coeff = v^Tc, since all numerical values of these parameters are 
known. Then the previous expansion looks like 



Pg{T) ~ ^A2(T2 + a;e//T)e-^iT _ _T [2?^ + 2ujeffT^ + '^e//^] 6""^^* 
16 



24, 



3 



^T^A^ 

— ^A^ 
27r2 



e//^ -r We// 

-J/3 



16 



2T^ + VSAT"^ + -A^T 
8 



(A36) 



The expansion (A36) clearly shows that the exponential suppression of any pressure at low temperature below Tc is 
determined by the corresponding asymptotic of the gluon mean number (4.6), namely 



Ng = ^—r^e-''^, Tc>T^O, (A37) 

CT — 1 

by replacing uj by h'Tc in each different case, as it is seen in the previous low-temperature expansion for the gluon 
pressure. For the scaled gluon pressure 3Pg{T)/T^ the expansion (A36) is especially useful, since it depends on the 
dimensionless variable {T/Tc) only, and it is shown in Fig. 1 below Tc. The expansion (A36) clearly shows that near 
to Tc the number of effective gluonic degrees of freedom and their magnitudes will be drastically increased, though 
the gluon pressure initially contains only five different massive and massless gluonic excitations. 



Appendix B: High-temperature expansion 



In order to evaluate a high-temperature expansion for the gluon pressure (Al), it is convenient to present the NP 
pressure (A2) as follows: 



Pnp{T) = A^T^ - ^A2Pi'(T) + ^TM{T), (Bl) 



since 



dw CO Ng{fi,(jj 



■ff 



Jo 



doj ujNg{(3,uj 



Jo 



dw w Ng{fi,(jj) 



(B2) 
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where 



/■oo /■oo 2 

duj w N,{P,oj) = da;^^^ = -n (B3) 



(B4) 



and the composition M{T) is aheady given by the relations (A4) and (A5). 

In the high-temperature Umit T — > oo (/3 = — > 0), the gluon mean number Ng{(3,oj) in the integral (B4) can 
be approximated by the corresponding series in powers of (/3a;), since the variable oj is restricted, i.e., 

NM^)^-^l-j={M-'[i-l{M + O{0% 13^0. (B5) 

Let us also note in advance that in what follows for our purpose it is sufficient to keep only the positive powers of T 
in the evaluation of the high-temperature expansion for the gluon pressure (Al), and hence for each term in Eq. (Al). 
Thus, for the asymptotic of the integral P{ (T) up to the leading order in powers of T, one obtains 



oo. (B6) 



In order to investigate the behavior of the composition M{T) (A4) at high temperature, it is convenient to decompose 
its integral P2{T), shown in Eqs. (A5), as follows: 

P2(r) = pi^)(r)-Pi')(r), (B7) 

where 

/•oo 4 2 

P^^\t) = r'"da;a;2ln(l-e-^"), (B8) 
Jo 



due to the relation (A30). Let us note in advance that we will not need the high-temperature asymptotic of the 
integral Pi^'(T). 

The integral P3(T) up to the leading order in powers of ^ — )• becomes 

P3(T) = ^ °" dw In (l - e'^'^'^ ~ ^ dw In ^w', /3 ^ 0, (B9) 

in accordance with the expansion (B5), since the variable io is restricted, and hence the variable lo' = \J liP- + a? as 
well, where a = v^A. The last integral can be exactly calculated and the high-temperature expansion for Pz{T) 
becomes 

P3(T) ^ ^-ulfj In - ^-u:lff + A^.;^// - V^^"" arctan , T ^ oo, (BIO) 

where ijj'^ff is given in Eq. (A12). 

The integral Pi{T) is convenient to decompose as the sum of two terms, namely 

P4(T) = Pi^\T) + Pf ^(T), (Bll) 
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where 



POO 

pf)(r) = / dw In (1 - e-'^") , 



(B12) 



Let us begin with the integral P^'^iT)^ which can be re-written as follows: 



poo 

) (T) = / dujujHn(l- e-/3-V(i+(«V4-^)) 



(B13) 



on account of the relation (3.6), namely w = \J ^jfi + (3/4)A2 = -^o;^ + (a/2)^. Since the variable oj is always big, 
then X = {a^ /Auj"^) <C 1, and thus we can expand 



^/{lTx) = l + ]^x + 0{x'^), x->0. 
Then the integral (B13) up to the leading order in powers of small /3 becomes 



(B14) 



dw a;2 In (l - e-/3-e-(-'3"/2)) , ^ ^ 0, 



(B15) 



where the argument of the exponent {xfiuj/2) = (a^/SwT) = 2 <C 1 in the T,u; — >■ 00 limit, so the integral (B15) can 
be present as follows: 



POO 

Pi^\T)r^ doj uj'^\n{l-e-^"[l- z + 0{z'^)]) , < 1, /3 ^ 0, 



(B16) 



or, equivalently. 



POO 

Pi\T) ~ / dww^ln 



P2(T) 



(1 _ e-/3-) 1 + 



eP" - 1 



dw In ( 1 



We// 



P2(T) + P^(r), 



(B17) 



as it follows from equations (A5) or (B7)-(B8). The argument of the logarithm in the second integral is again always 
small {zje^^ — 1) <C 1 in the T, a; — )• 00 limit, and thus we can expand it and obtain in the leading order 



- 1 



(B18) 



The last integral is nothing but Pi{T) defined in Eq. (B2), so that combining (B2) and (B6), one obtains 



P^(r)~-/3Pi(r)~-A2T--AV//, T^oo, 



(B19) 



and then the high-temperature expansion for pj^^ (T) becomes 



P«(r)~P2(T) + ^A2r-^AV//, T^oo. 



(B20) 
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("2) 

The integral P4 ' (T) in the leading order in powers of /3 — )• becomes 

pP{T) = dujLj'^ In (1 - e-f^^) ~ / dw i^^^^ ^ ^ (B21) 

Jo Jo 



in accordance with the expansion (B5), since the variable w is restricted, and hence the variable w = ^/ui"^ + 

as well, where again a = \/3A. Prom the relations (3.4) and (3.5) it follows that w' ^ a) by A ^ (1/2) A, so making 

this replacement in the expansion (BIO), one automatically obtains the high-temperature expansion for the integral 

Pi^\T) as follows: 

Pi'' - l-hf m (1^) - l-hf + \^'-eff - f A3 arctan (^) , ^ ^ 00, (B22) 

where uj^jj = Loljf + (3/4) A^. 

The high-temperature expansion for the composition (A4), on account of the relations (B11)-(B12) and the previous 
expansions (BIO), (B20) and (B22) and after doing some algebra, becomes 

l^TM(T) ^ l|AWT-A^T^+3^.3^,Tln(<£) 

+ ^A^Tarctan ( - ^A^T arctan ( ^\ , T ^ 00. 

Substituting this expansion into the Eq. (Bl), and on account of the expansion (B6), we obtain 



(B23) 



+ ^A^T arctan ( - l^A^T arctan ( , T ^ 00. 



(B24) 



One concludes that the exact cancelation of the A^T^ term occurs within the NP pressure itself. Thus Pnp{T) ~ T up 
to the leading order in the T ^ 00 limit. At the same time, at high temperatures the exact cancelation of the Psb{T) 
term occurs within the composition {\Q / tt^)T Mi{T) = {IQ / ■k'^)T[P2{T) ~ Pi{T)], which enters the composition (B.23). 
To show this explicitly, let us substitute into the former composition the relation (B7), on account of the relations 
(B8), and the relation (Bll), on account of the expansion (B20), and doing some algebra, one obtains 

i^TMi(T) ~ -2Psb{T) + 2Psb{T) - A^T^ + ^A^uj^ffT - ^^TPf ^(T) 

^ _A2t2 + AAV//T-^rpf' (T), T^oo, (B25) 

TT TT 

from which the above-mentioned exact cancelation explicitly follows. The exact cancelation of the ^2^' (^) terms and 

(2) 

the expansion (B22) for the P^ (T) term arc not shown, for simplicity. 

Let us now consider Eq. (6.4), which is convenient to decompose as follows: 



PpAT) = ^«.A2 / dco co' -^^j^ = PliT) - P|(r), (B26) 



where 
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1 1 



Ztt^ Jq uj c 1 



(B27) 



and let us remind that w = i/w^ + (3/4) A2 = -y/co^ + (a/2)2. 

In the integral Pl{T) it is convenient to introduce a new dimensionless variable x = j3uj = + {a/2)^. After 

doing some algebra, it becomes 



Pl{T) = ^asA' / doj oj' = ,r2^sA^T^ / dx 



v/:c2 _ (/3a/2)2 



- 1 



The last integral when ^ — )• can be approximated up to the leading order as follows: 



dx 



{Pa/2) 



- 1 



- 1 6 



then for the integral (T) up to the leading order in powers of T we get 

PUT) ^^asA'^T^, T^oo. 



(B28) 



(B29) 



(B30) 



In the integral PUT) the variable u is restricted, and hence to as well. So up to the leading order in the T 
00 = 0) limit this integral can be approximated as 



P^(T) 



duj — - 



w eP'^ - 1 27r2 



AyM ^2 

ra^A^T / dco T ^00, 



(B31) 



in accordance with the expansion (B5). The last integral can easily be integrated and thus the high-temperature 
expansion for the P^iT) term looks like 



PUT) ^ ^ccsA^T 



V3 ^ (2Kym\ 
— A arctan — = — 



AvM ^ A arctan 



2 V V3 A ; 

Summing it with the expansion (B30), for the integral (B26) one obtains 



, T 00. 



(B32) 



Ppt{T) ~ ^asA^ 



-T^-T A.M-^Aarctan(-^ 



T 00, 



(B33) 



which is nothing but the high-temperature expansion for the a^-dependent PT part of the gluon pressure. 

The high-temperature expansion of the gluon pressure is to be obtained by summing up the expansions (B24) and 
(B33), so it is 



P,{T) = [P^p(r) + P^^(r)]~^AV//r+^a;,^//rin(g^) 
+ ^A^T arctan ( - l^A^T arctan ( 



9 a2 



V ySA y 7r2 
-T"^ -T\ Aym - -V A arctan 



\V3AJ 

2AyM \ 

V3A J 



T 00. 



(B34) 
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Let us emphasize that the high-temperature expansions for ah three pressures (B24), (B33) and (B34) non-analyticahy 
depend on the mass gap in terms ^ A'^T ^ (A^)^'^/^)r, but it is not an expansion parameter hke Us- From asymptotics 
(B33) and (B34) it fohows that Pprp{T), and hence Pg{T), behaves hke in the leading order, while remembering 
Pnp{T) ~ T, see expansion (B24). It is also interesting to note that the effective massive gluonic "excitations" 

Wg^j = \J^'iff + 3A2 and We// = \J^1ff + (3/4)A2 are logarithmically suppressed at high temperatures, while there 

is no dependence on the effective massive gluonic " excitation" Co^f f = \/ ^^-^ + (3/4) A^ at all. 
In a more compact form the previous expansion looks like 



Pg(T) - Baa^A^T^ + [BgA^ + Af3]T, T ^ oo, (B35) 

where denotes the terms of the dimensions of the GeV'^, which depend analytically on the mass gap A^. The 
explicit expressions for it and for both constants B2 and -B3 can be easily restored from the expansion (B34), if 
necessary. 
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